6 - Introduction to the Finite Element Method - Tutorial [ID:52629]
50 von 367 angezeigt

Yeah, so welcome to our today's tutorial.

Today we will consider exercise four. I don't know what's wrong here.

Hope it remains as it is. So exercise four. Here we shall a little bit practice the integration

with different quadrature rules and to this end we have given three different functions.

They are all polynomials and we should do a numerical integration using either Newton

codes quadrature or Gauss Legendre quadrature. So let me just copy this.

So here we go. So in the very first part here we have to integrate the function f1 using Newton

Gauss quadrature, Newton codes sorry, quadrature with two quadrature points which is the trapezoidal

rule and let me furthermore open the lecture notes because there we have this nice overview

of the weights and gratitude points. So this is here. So this we can keep in mind here.

Okay. So the first function here. This is a linear function so you can also say a polynomial

first order and first let us compute the analytical solution because this is also

asked for. So the analytical solution here. It's just a standard integration. The limits

are given here. It's the interval from minus one to one. So what do we get after integration here?

Xi squared plus I'm not happy with just Xi squared. If you take the derivative

of Xi squared what do you get? Yes that's it. That's what I wanted to know. Slight different

result could follow here. Okay so what do we get? One divided by two? No.

So it's actually

so now it should yes it's two.

When we discuss the accuracy of Newton-Coltz gratitude we distinguish between even and odd

polynomials. Do you see what could be the reason for that? Here by the analytical treatment. What

happens to this term where we have after integration? Exactly it cancels out because if you take

one and minus one in general to an even power then it just cancels out in this interval from

minus one to one. So and the even exponent here actually results from an odd exponent before

integration. So actually the odd exponent does not play a role. Okay so this is the analytical

solution and now we can integrate with Newton-Coltz and we choose two integration points.

Would you please look up which degree of polynomial is possible to be integrated

exactly by the choice here. Just have a look to your notes. The morning we discussed Gauss

quadrature. This was this formula two times the number of quadrature points minus one. This holds

for Gauss. N minus one. Why N minus one? Yeah maybe let's have a look here. This was in notes of last

time. So

just as a side remark here.

So if we have an even polynomial then yeah it's two minus one. If we have an odd polynomial it's

two but two is not odd. So in fact we can integrate by these two

quadrature points a linear function exactly. So and let us try whether this works. It should work of course.

So the approximation is that we

evaluate the first node and then we evaluate the function at the second node, sorry second quadrature point.

So which is in our case

to be evaluated at

here at minus one and at one and the associated weights are one in both cases. So we have to write here

the function one by two evaluated at minus one multiplied by one plus

so I put it here in a very formal way of course one can abbreviate that and what is the result here?

What do we get? Two. So as expected. Okay

so NEWTON CODE with two integration points

what do we get? 2. so as expected. okay so Newton codes with two integration

points

so based on this finding we can have a look at the next function which is here

a cubic function yeah and for sure with the two so here we are 4.2 with the two

integration points we will for sure not be able to integrate it exactly. how many

quadrature points do we need to integrate a cubic polynomial exactly?

Zugänglich über

Offener Zugang

Dauer

01:31:19 Min

Aufnahmedatum

2024-05-07

Hochgeladen am

2024-05-09 02:39:33

Sprache

en-US

Einbetten
Wordpress FAU Plugin
iFrame
Teilen